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ABSTRACT 

We present observations of the nuclear star cluster in the nearby edge-on spiral galaxy NGC 4244 
using the Gemini Near-Infrared Integral Field Spectrograph (NIFS) with laser guide star adaptive 
optics. From a previous study of edge-on galaxies, this nuclear star cluster was found to be one 
of a sample of clusters that appear flattened along the plane of their host galaxies disks. Such 
clusters show evidence for multiple morphological components, with younger/bluer disk components 
and older/redder spheroidal components. Our new observations of NGC 4244 show clear rotation of 
30 krns" 1 within the central 10 pc (0.5") of the cluster. The central velocity dispersion is found to be 
28±2 kms -1 . The multiple stellar populations inferred from the optical colors and spectra are seen 
as variations in the CO line strength in the NIFS spectra. The rotation is clearly detected even in the 
older, more spheroidal stellar component. We discuss evidence for similar structures and kinematics 
in the nuclear star clusters of other galaxies including M33 and the Milky Way. Our observations 
support two possible formation mechanisms: (1) episodic accretion of gas from the disk directly onto 
the nuclear star cluster, or (2) episodic accretion of young star clusters formed in the central part of 
the galaxy due to dynamical friction. 

Subject headings: galaxiesmuclei - galaxiesrstar clusters - galaxies: spiral - galaxies: kinematics and 
dynamics - galaxies:formation - galaxies: individual (NGC 4244) 



1. INTRODUCTION 

Prominent nuclear star cluster s (NSCs) are foun d in 
~75% of late- type s piral galaxies dBoker et al.ll2002f ) and 
early-type galaxies (jCote et aT1l2006| ). These NSCs are 
intrinsically very luminous, with typical Mr ~ —12, and 
sizes similar to globular clusters {r e t f ~ 5pc;|Boker et al.| 
2004) . Recently, multiple studies ilFerrarese et al.1 2006 



from the observations of several distinct populations in 
the color-magnitude di agram of the nucleus of the Sgr 
dwar f spheroidal (M54; iMonaco et~aTll2005l: ISiegel et all 
[200% 



Wehner fc Harris! 12009 : iRossa et all I2Q06T ) have shown 
that NSC masses are coupled to the mass of their host 
galaxy, following a similar M-a relation as massive black 
holes (MBHs) and extending this relation to lower mass 
galaxies. This result suggests that the formation of both 
types of objects are linked to the formation of their host 
galaxies. Some galaxies host both NSCs and MBHs, but 
the relati onship between t hese two types of objects is 
unknown (Set h et al]|2008f h Unlike black holes, the as- 
sembly history of NSCs can, in principle, be inferred from 
studies of their stellar populations and kinematics. Thus 
NSCs provide important tools for studying the history of 
galaxy centers and their connection to the galaxies as a 
whole. 

In spiral galaxies, integrated optical spectra of NSCs 
show tha t they nearly always have multiple stel lar pop- 
ulations (|Walcher et all 120061 : IRossa et al.ll2006[ K More 
specifically, NSCs typically have a young stellar popu- 
lation (<100 Myr) but are dominated in mass by stars 
older than >1 Gyr. Further support to this picture comes 
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In ISeth et all (|2006l ) (Paper I) , we found from a study 
of edge-on galaxies that in some NSCs these multiple 
stellar populations appear to be associated with different 
morphological components. These clusters are flattened 
along the major-axis of the galaxy and have color maps 
and morphologies that suggest they consist of young 
disks and older more spheroidal components. 

We present here follow-up observations of the NSC in 
NGC 4244, the nearest galaxy in the P aper I sample 
(D=4.3 Mpc, 1"=21 pc: ISethet all 12005?) . In Paper I, 
we presented integrated optical spectra of this cluster 
that indicate the presence of multiple stellar populations. 
Although the detailed star formation history is not well- 
constrained, populations of 50-100 Myr combined with 
older populations are required to provide a good fit to 
the spectra. Redshifted emission identified with an HII 
region offset by 19 pc along the major axis of the cluster 
gives a lower-limit to the mass of ~2.5xl0 6 M Q . This 
mass measurement is consistent with the masses derived 
from the stellar population fitting of the spectrum as 
well. 

In §2 we describe the observations and data reduction. 
Our results on the kinematics, morphology and stellar 
populations of the cluster are presented in §3. We discuss 
the results in the context of other galaxies and formation 
models in §4. 

2. OBSERVATIONS & DATA REDUCTION 

K-band spectra of the NGC 4244 nuclear star clus- 
ter were obtained using the Near-Infrared Integral Field 
Spectrograph (NIFS) on Gemini North, an image-slicing 
integral field unit spectrograph. The spectrograph was 
used in conjunction with the ALTAIR Laser Guide Star 
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system. The cluster itself was used as a tip/tilt star, 
but is too faint (r=16.18) to get the highest-Strehl per- 
formance. The cluster was observed in K-band on the 
nights of February 9 and April 29 2007. In total 5x900 
second exposures were obtained on-source and 4x900 sec 
exposures off-source (sky). The airmass of the obser- 
vations ranged between 1.11 and 1.42, and an A0 star, 
HIP50459 was observed on both nights at similar air- 
mass to be used as a telluric calibrator. The exposures 
all had similar image quality, with the PSF FWHM vary- 
ing by <30% between the exposures. We note that ri- 
band spectral observations of the NGC 4244 NSC were 
obtained as well, but due to a number of observational 
problems and a significantly lower resolution, these data 
are not discussed here. 

The data were reduced using the Gemini v. 1.9 IRAF 
package, utilizing pipelines based on the NIFSEXAM- 
PLES scripts. To obtain a clean telluric absorption spec- 
tra from the HIP50459 observations, the spectra were 
processed as described for the science images below, after 
which 1-D spectra were extracted and the Br7 absorption 
line was fit and removed. 

For the science exposures, an off-source exposure was 
subtracted from each on-source exposure; for the Febru- 
ary 9 data only one off-source exposure was available 
and was thus subtracted from both images. The images 
were then sliced up and rectified based on Ronchi mask 
and arc lamp images using NFTRANSFORM and spa- 
tially rebinned into image cubes using NIFCUBE. In this 
process, the original 0.043x0.1" pixels get rebinned into 
0.05" pixels which are block replicated in the long (slice) 
direction. Each data cube was then corrected for telluric 
absorption using the NFTELLURIC program, which al- 
lows for interactive correction of the small differences in 
the wavelength calibration and airmass between the ob- 
servations of the calibrator and each exposure. The fi- 
nal data cubes from each individual on-source exposure 
were combined using GEMCOMBINE correcting for spa- 
tial and spectral offsets. The resulting data cube has a 
-3"x3" field of view with 0.05" spatial pixels and 2.132 
A spectral pixels from 2.00-2.44 /im. Based on measure- 
ments of the sky lines, the spectral resolution is ~4.lA 
(53 kms -1 ) FWHM. The spatial resolution in the com- 
bined image was measured using four compact sources 
after subtracting the NSC model (§3.1) and found to be 
0.23" FWHM. These four compact sources are shown as 
crosses in Fig. [1] three of the four are also detected as 
point-like sources in the higher resolution (0.08" FWHM) 
HST images. The spatial resolution of our data cube 
is roughly ^2x lower than that achieved when using a 
brighter tip/tilt star (e.g. for our telluric calibrators). 

A spectrum from the inner 0.1" of the cluster is shown 
in Fig. [TJ The plotted spectrum is the normalized flux 
ratio of the galaxy spectrum to the Br7-corrected A0V 
telluric calibrator spectrum. All the spectral analy- 
sis presented below is done using flux ratioed spectra. 
In order to achieve a reasonable S/N for analysis of 
the outer parts of the cluster, the data were spatially 
bin ned using the Voronoi tes selation method described 
in (|Cappellari fc Copinll200l . The S/N of each pixel 
was determined empirically using a relatively featureless 
portion of the spectra between 22140 and 22866A. Indi- 
vidual pixels near the center of the image have S/N>50 
per pixel. We binned pixels together in order to obtain 
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Fig. 1. — Top - The spectrum within 0.1" of the nuclear star 
cluster center. The plotted spectrum is the normalized flux ratio of 
the galaxy spectrum to the Br7-corrected A0V telluric calibrator 
spectrum. The CO bandheads dominate the red part of the spec- 
trum while Mg I (21080,228301), Br7 (21676A), Na I (22075A), 
and Ca I (22644-77A) lines are also seen. Bottom - K-band image 
made from collapsing the final data cube along the spectral direc- 
tion. Contours are at 16.5,15.5,14.5, 13.5, and 12.5 mag/arcsec 2 . 
Crosses show the four compact sources used to determine the spa- 
tial resolution of the data. The compass in the upper left shows 
the rotation of the field (PA=147°); the x-axis runs roughly along 
the major axis of the galaxy. 

a S/N~25 in each bin, so near the center, the 'bins' con- 
sisted of a single pixel. We found that the S/N did not 
increase as expected with binning (at least in part due 
to spectral rcbinning during the initial reduction), thus 
the outer bins have S/N of <10. 

We also reconstructed a calibrated K-band image of the 
nuclear star cluster for the morphological analysis pre- 
sented in §3.1. To do this, we first multiplied the spectra 
by a 9500K black body to correct for the division by the 
A0 telluric calibrator, and then multiplied this spectrum 
by the 2MASS K s filter response curve before collapsing 
the data cube to an image. We then calibrated the zero- 
point of this image by comparing the flux from the best 
fitting model (described below) to the magnitude of the 
source in the 2 MASS Point Source Catalog (iiT s =13.55; 
iCutri et al.ll2003f ). Based on the variation between mod- 
els and the uncertainty in the sky-level on our image, this 
zeropoint should be good to ~10%. The image and cali- 
brated contours are shown in the bottom panel of Fig. [TJ 

3. RESULTS 

In this section, we present results from our Gem- 
ini/NIFS K-band observations. We first discuss morpho- 
logical fits to our high resolution K-band image (§3.1). 
We then present the kinematics of the cluster in §3.2 
and its stellar populations in §3.3. The most significant 
finding is that the NGC 4244 NSC is strongly rotating. 
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TABLE 1 
Results of Morphological Fits 

NSC Disk Component Fits 





M> [Lq/p^] 


h x [pc] 


20 [pc] 


K-band disk+Sersic 1.4xlO b 3.39 
K-band disk+King 1.9xl0 5 2.72 
F814W disk+King 9.7xl0 4 2.73 

NSC Spheroidal Component Fits 


1.19 
1.38 

1.42 




M> [L©/pc^] 


r eff M 


(b/a) 


K-band disk+King 
F814W disk+King 


4.5xl0 4 
5.3xl0 4 


13.6 
5.70 


0.79 
0.73 




Peff [ L ©/P c ^] 


r eff [pc] 


(b/a) n 


K-band disk+Sersic 


8.7xlO a 


10.9 


0.81 1.68 



Note. — All quantities derived assuming a distance of 4.3 Mpc. 

Comparison of the kinematics to the stellar populations 
and morphological fits shows that this rotation is occur- 
ing even off the major-axis of the cluster where older 
stellar populations appear to dominate. 

3.1. Morphological Results 

We briefly describe our analysis of the morphology of 
the NGC 4244 NSC using the high resolution K-band im- 
age of the cluster created from our Gemini/NIFS data. 
This information is used in §3.3 to argue that the entire 
cluster is rotati ng, and will be used as input for dynam- 
ical modeling flde Lorenzi et al.l l2008f) using NMAGIC 
code (|de Lorenzi et al.ll2007D . 

As we found for the optical data (see Paper I, Fig. 4), 
a single component elliptical King (or Sersic) profile fit 
to the K-band image leaves behind significant residuals 
in the plane of the galaxy. Combined with the opti- 
cal color difference (Fig. [3] top panel) , this motivates us 
to fit the data using both a spheroidal and disk com- 
ponent. All models were convolved by a Gaussian PSF 
with FWHM of 0.23" during the fitting process. We tried 
both disk+King (as used in Paper I) and a disk+Sersic 
profile. Both profiles provide fits with reduced x 2 values 
~4x smaller than the spheroidal component alone. The 
values of the derived parameters for these fits, and the 
disk+King F814W filter fit from Paper I, are shown in 
Table [T] The parameters for the disk fits are the cen- 
tral surface brightnesss /j,q, the scale length h Xl and scale 
height zq. For the spheroidal component fits, the pa- 
rameters are the central surface brightness [1$ (King) or 
effective surface brightness fJ-eff (Sersic), effective radius 
r e ff, axial ratio (b/a), and Sersic index n. The derived 
K-band disk sizes are very similar to those derived at op- 
tical wavelengths from the HST data despite the lower 
resolution of this data. The spheroid geometry is signifi- 
cantly more extended in the K-band, but its size is quite 
sensitive to the fitted value of the background. 

Based on disk fits t o the full galaxy from 2MASS pre- 
sented in ISeth et all (|2005h , the central luminosity sur- 
face density of the NSC is ~250x greater than that of 
the disk, while the luminosity surface density at the edge 
of the NIFS field is comparable to the level predicted by 
the galaxy disk model. Deprojecting the disk compo- 
nents, the luminosity volume density at the center of the 
NSC is ^10 6 x that of the underlying galaxy disk. The 
local volume density of the cluster thus dominates the 
mass in the central regions of the galaxy out to the edge 
of our data. 



3.2. Kinematic Results 

We first describe the two ingredients used to derive 
the kinematics of the NGC 4244 nuclear star cluster: 
(1) spectral templates, and (2) a method to determine 
the line-of-sight velocity distribution (LOSVD) from the 
binned spectra based on the templates. 

Spectral templates were obtained from two sources. 
Three templa t e sta rs were observed with NIFS by 
iDavidge et al.l (|2008| ) in a configuration identical to ours. 
Spectra for these stars were kindly provided by the au- 
thor. Their spectral types: K1III, K5III, and M0III 
are representative of the spectral types of the red gi- 
ant branch (RGB) and asymptotic giant branch (AGB) 
stars most likely to dominate our spectrum. We also 
use a larger library of spectral templates that has been 
observed by C. Winge using the somewhat higher reso- 
lution (3.35 A FWHM) GNIRS spectrograph on Gemini 
South 7 . This library contains stars from spectral type 
F7-M3 and includes dwarf, giant and some supergiant 
stars. 

To determine the LOSVD of each spectrum we 
used the penalized pi x el fit ting (PPXF) method of 
ICappellari fc Emselleml (|2004D . The program finds the 
best matching set of templates and convolves them to 
find the best fitting LOSVD, which is parametrized us- 
ing a Gaus s-Hermite series. The non-gau ssian terms (ft.3 
and M, see Ivan der Marel k. FrTanxll 19931 ) in the LOSVD 
are penalized, thus the solution typically converges to 
a Gaussian LOSVD at low S/N. The method also com- 
pensates for any low-frequency differences between the 
template and galaxy spectral shape using Legendre poly- 
nomials. We fit the area of the spectrum from 22900A to 
24000A, which is dominated by strong CO lines. We cal- 
culated errors on the LOSVD parameters using Monte 
Carlo tests that consisted of adding noise to the spec- 
trum based on the measured S/N and remeasuring the 
LOSVD 250 times without penalizing non-Gaussian dis- 
tributions. The typical standard deviation on the veloc- 
ity and dispersion for pixels with S/N~25 is ^3 kms -1 . 
When using GNIRS templates, we convolve the template 
spectra by a Gaussian with a — 13 kms -1 to correct for 
the slightly higher resolution of the GNIRS data. 

We show the radial velocity map derived using PPXF 
on each spectral bin in the top panel of Figure [2] We 
show only data with reliable measurements by clip- 
ping bins with S/N<10 or errors on the velocity above 
15 kms -1 . The radial velocity map shows strong rotation 
in the NGC 4^44 NSC The rotation velocities reach an 
amplitude of 30 kms -1 ^0.4" from the nucleus. Assum- 
ing a Keplerian disk, this rotation suggests a mass of 
1.7xl0 6 solar masses enclosed within the central 8 pc. 
Beyond this radius, there is some suggestion the velocity 
may drop, however due to the low S/N at large radii this 
result is not robust. 

Based on the morphological results presented in the 
previous section, the NSC dominates the mass in the 
center of NGC 4244, and thus is likely to be kinemati- 
cally d i stinct from the rotation of the overall galaxy disk. 
lOllingi (|1996f ) has presented a study of the HI rotation 
curve of NGC 4244 with resolution as high as 10"x20". 
The HI gas rotates in the same direction as the NGC 4244 

7 http://www.geniini.edu/sciops/instruments/nir/spectemp/index.htnil 
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NSC. The innermost point in their derived rotation curve 
is at 24" where the galaxy rotation is 35 kms -1 , increas- 
ing to 60 kms -1 at 47". Given that we observe rotation 
of ~30 km s -1 within the central half-arcsecond, this sug- 
gests that the rotation of the NSC is distinct from the 
gas component of the galaxy disk. It would be interest- 
ing to compare the NSC rotation to the stellar rotation 
in the galaxy disk, however, this would require relatively 
deep long-slit spectroscopy (e.g. lYoachim fc Dalcantonl 
l2008f) since the surface brightness of the galaxy beyond 
the nuclear star cluster is quite low. 

The line-of-sight velocity dispersion is a somewhat 
more difficult quantity to measure accurately due to the 
low dispersion of the cluster as a whole. Our spectra 
have 2.1A/28 kms -1 pixels and 4.1A/53 kms -1 reso- 
lution (FWHM), while the dispersions of NSCs in late- 
type galaxies are typically ~25 kms -1 (jWalcher et al.1 
2005). Nonetheless, we are able to make an accurate 
measurement of the central dispersion of the cluster, 
a = 28 ± 2 kms -1 . Convolution of the templates with 
this dispersion substantially (by a factor of 2) improves 
the x 2 °f their fit to the central spectrum of the cluster, 
giving a reduced \ 2 0I 0.9 and 1.1 using the GNIRS and 
NIFS templates respectively. The bottom panel of Fig. [5] 
shows a map after clipping data with S/N<15 or errors 
in the dispersion >10 kms -1 . The velocity dispersion 
drops along the major axis to to values of 10-20 kms -1 
at precisely the area where rotation is seen most strongly 
indicating that we are indeed seeing a relatively cold ro- 
tating disk. This dispersion is close to the limit of what 
we can reliably recover from our relatively low resolu- 
tion observations, thus these values may represent upper 
limits. It is unclear whether the increased dispersion we 
see in the center is a result of unresolved rotation or a 
genuinely hot component. 

The use of the templates from NIFS and GNIRS give 
very similar results for both the velocities and disper- 
sions. The GNIRS templates selected for fitting the spec- 
tra are primarily at spectral types close to those (K-M 
giants) represented by the 3 NIFS spectra. For almost 
all the spectra, our measured non-Gaussian /i3 and M 
components of the LOSVD were consistent with zero. 
However, due to our low S/N and resolution, we would 
not ex pect to be able to detect any non-Gaussian signa- 
tures (ICappellari fc Emsellemll2004D . 

3.3. Stellar Populations and Connection to Kinematics 

From the optical data described in Paper I, the 
NGC 4244 NSC has two distinct stellar populations - a 
bluer population dominates in the midplane with a red- 
der population dominating above and below the plane. 
This is clearly seen in the HST based color-map shown 
in the top panel of Fig. [3] We would expect the NIFS 
data to resolve these different populations despite the 
2x lower resolution. If we make a map of the CO line 
strength (bottom panel of Fig. [3]) , the bluer areas of the 
HST image clearly have stronger CO absorption than the 
redder areas above the midplane. The CO line strength 
is determined by comparing the average depth of the line 
at rest wavelengths of 22957±52A to the continuum mea- 
sured just to the blue of the CO bandhead ( 12 CO (2, 0); 
iKleinmann fc Halj [l986). We also measured these values 
for the GNIRS template stars for comparison. The CO 
line strength values for the NSC range from ^0.10 above 
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Fig. 2. — Top - Color indicates the measured radial velocity 
for the Voronoi binned data using templates observed with NIFS. 
Gray areas indicate spectra with S/N < 10 and/or errors > 15 
kms -1 . Rotation of ~30 kms -1 is clearly visible along the ma- 
jor axis. Contours show the K-band image at 15.5,14.5,13.5, and 
12.5 mag/arcsec 2 . The black bar indicates 10 pc (0.47"). Bot- 
tom - Velocity dispersion measurements with S/N >15 and errors 
<10 kms -1 . The dispersion drops away from the center indicating 
a relatively cold disk population. 



and below the major axis typical of late G-type giant 
stars (T e ff ~ 4800 K), to ~0.25 along the major axis, 
typical of K3-K6 giants (T eff ~ 4100 K). We caution 
that the CO depth depends on the gravity of the stars as 
well as their effective temperature. We note that while 
Mg, Ca, and B1-7 absorption are clearly detected in the 
highest S/N spectra, these features have lower S/N than 
the CO feature and do not show any obvious variations 
across the NSC. No emission features are detected. 

To better understand the relation between the optical 
colors and CO line strengths seen within the cluster, in 
Fig. 0] we model these quantities using Padova isochrones 
which include an updated treatment of the asymptotic gi- 
ant branch (AGB) (Marigo et al. 2008). We sampled 10 7 
Mq of stars from a iKroupal (|2001[ ) initial mass function 
and then obtained effective temperature, gravities, and 
the F606W, F814W, and K-band magnitudes for each 
star from the isochrones at ages from 4 Myr to 14 Gyr. 
To obtain estimat es for the CO line strength, we used 
the relations from Dov on et al.l (|1994f ). who fit CO line 
strength as a function of effective temperature for dwarfs, 
giants, and supergiants. The total CO line strength at 
each age was calculated as the mean CO line strength 
weighted by the K-band lumin o sity o f each star. We note 
that the broader iDovon et~a l. (1994) CO index has val- 
ues that are 20-50% lower than the one we measure here, 
nonetheless, the trends of CO depth with age should be 
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Fig. 3.— 2"x2" images of top - the F606W-F814W color map 
from the HST/ACS data described in Paper I. The image has been 
rotated to match the orientation of the NIFS data. The contours 
indicate the F814W surface brightness and were chosen to roughly 
match the K-band contours in [2] and in the bottom panel. The 
black bar indicates 10 pc (0.47"). Bottom - the CO line strength 
with larger values indicating deeper CO lines. The contours show 
K-band contours as inl2l 

robust. We also caution that the models remain very 
uncertain, especially given that the CO strength can be 
dominated by the coolest, thermall y pulsing AGB stars , 
which are very difficult to model (|Marigo et al.l [2008) . 
The CO line strength shows peaks at two ages; first at 
10-20 Myr due to the presence of bright red supergiant 
stars, and second at ~300 Myr from AGB stars. The 
CO evolution of these mo dels is qual i tative ly similar to 
that seen in the models of iMarastonl (|2005f ). which also 
include contributions from thermally pulsing AGB stars, 
but differs from earlier models that don't include this 
phase and thus are l acking the second peak in CO depth 
(e.g. iLeitherer et al.lll999f) . 

We can use these models to qualitatively analyze the 
variations in stellar population seen across the cluster in 
Figs. [3] and [5l The extinction is not well known, however 
spectroscopic fits suggest an Ay ^0.5 corresponding to 
a reddening of ~0.2 magnitudes in the F606W-F814W 
color. For the youngest stars, we expect a metallicity 
of Z^O. 008 based on the 100 krns" 1 circular velocity of 
NGC 4244 (|Garnettll2002l ). However, the combination of 
the red color and low CO values above and below the 
midplane of the cluster suggest an old (>1 Gyr), more 
metal-poor population. Along the major axis there is a 
slight reddening in the optical color along with an in- 
crease in the CO line strength. One possible explanation 
for this is an age gradient within the disk; for instance 
if the center had an age of ~100 Myr years then the 
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Fig. 4. — Mode ls of the optical colo r CO line strength using 
Padova isochrones IMarieo et all (2008) at Z=0.004, 0.008, and 
0.019. The to p panel shows the CO line strength derived using 
relations from Dovon et al. (1994), while the bottom panel shows 
the F606W-F814W optical color. 

redder colors and deeper CO could result from slightly 
older or younger populations at 5-10 pc. We caution 
that the stellar populations at any point in the cluster 
are clearly mixed and thus information determined from 
comparing the observations to single-stellar population 
models should be viewed with caution. With the ma- 
turing of NIR stellar models and spectral librar ies (e.g. 
iLanconet aIll2007t[Mlrmol-Queralt6 et al.ll2008l ). future 
population synthesis studies may be able to provide more 
robust results on the resolved star formation history from 
this type of data. Lastly, we note that it is possible for 
thermal emission from warm dust to reduce the strength 
of the CO absorption. However, the low star formation 
activity in the NSC indicated by the lack of NIR emis- 
sion lines suggests that dust emission probably does not 
play a large role in this cluster. 

We now turn to the relation between the stellar pop- 
ulations and the rotation. Figure [5] shows the derived 
velocity, CO line strength, optical color, and disk frac- 
tion of the NGC 4244 NSC as function of position in 
slices along, above, and below the major axis. The slices 
are each 0.25" in height, and the plotted lines are the 
result of averaging the maps (i.e. Fig. [5] and [3]) along the 
vertical direction in each slice at every pixel along the x- 
direction. The velocity and CO measurements were made 
using the Voronoi binned data; the symbols in these pan- 
els show the light- weighted positions of the bins coded by 
which slice they fell in. 

The top panel of Figure [5] shows that the rotation is 
not confined to the midplane, but is present ~0.25" above 
and below the midplane as well. The bottom three panels 
show that the slices above and below the plane are clearly 
sampling a different component/stellar population from 
the midplane. The CO lines are less deep away from the 
plane, while the optical color is redder. The disk fraction 
in the bottom panel is based on the results of our mor- 
phological descomposition from §3.1. We calculated the 
fraction of disk light along our slices from the best-fitting 
Sersic+Disk model after convolving each component with 
the PSF. For the area above and below the major axis, 
only 10-20% of the light in our NIFS observations appear 
to be coming from the disk. Taken together, these plots 
show that the whole cluster is rotating, not just the bluer, 
younger disk component. This finding places strong con- 
straints on the formation of this NSC as detailed in the 



Seth et al. 



Major Axis Offset [arcsec] 
-0.4 -0.2 0.0 0.2 0.4 




-15 -10 -5 5 10 15 

Major Axis Offset [pc] 

Fig. 5. — Profiles of (from top-to-bottom) the relative radial ve- 
locity, CO line strength, optical (F606W-F814W) color, and frac- 
tion of light expected from the disk based on the morphological 
models. Fach quantity is plotted for the midplane (black solid line 
and symbols) and for slices above and below (red dashed line and 
symbols) the major axis. Lines show vertically binned averages in 
each slice, while the symbols in the top two panels show the values 
of each bin that fell within our slices. 

following section. 

4. DISCUSSION & CONCLUSIONS 

In NGC 4244, we find a flattened nuclear star clus- 
ter with multiple components, a younger disk structure 
and an older spheroidal structure (Paper I). From our 
NIFS data, we have shown that the cluster is rotating at 
±30 kms -1 within the central 10 pc of the cluster. This 
rotation is not confined just to the younger disk compo- 
nent, but also appears to be present in the older, more 
spheroidal component. 

There are two basic models for NSC formation: (1) 
globular clusters which collect at the center due to dy - 
namical friction (|Tremaine et al.|[l975t lLotz et al.ll2001f ). 
or (2) episodic accretion of gas onto the center fol- 
lowed by in situ star formation (Milosavlicvic 2004; Bckki 
120071 ). A third alternative (3), combining these two sce- 
narios, has recently been proposed by Milosavljevic & 
Agarwal (2008, private communication) in which NSCs 
form from young star clusters in the inner part of the 
galaxy disk which are massive enough to accrete onto 
the nucleus. These mechanisms are not exclusive; NSCs 
could form from a combination of these effects, and differ- 
ent types of galaxies may have different dominant mecha- 
nisms. However, our observations for NGC 4244 strongly 
conflict with the first scenario. If a significant fraction of 
the NSC were formed from globular clusters distributed 
in a spheroidal distribution no rotation would be ex- 
pected. Thus the primary formation mechanism must 



be through episodic accretion of material from the disk 
of the galaxy, either gas or young star clusters. Distin- 
guishing between gas and star cluster accretion may be 
possible through comparison of our kinematic data with 
simulations. Direct evidence of gas accretion into the 
central ~10 parsecs is seen in a co uple nearby nuclear star 
clust ers using CO observations (|Schinnerer et all 120031 . 
[200(1 . 

There is substantial evidence that the structure we see 
in NGC 4244 is not unique, and in fact may be typical for 
NSCs in spiral galaxies. In Paper I, 5/9 NSCs in edge-on 
spirals were significantly flattened along the major axis 
of the galaxy disk, with three showing clear composite 
morphologies. NGC 4244 is the nearest of these galax- 
ies and thus presents the clearest view of these struc- 
tures. There is also evidence for s i milar structures in 
earlier type galaxies; iBalcells et al.l (|2007l) find a num- 
ber of disky "nuclear extended components" with scale 
lengths as small as 5 pc in a sample of early-type spirals 
and elliptical gala xies. Previous studies of nuclei in large 
samples of spirals (Baker et al. 2002; Carollo et al. 2002) 
have focused on face-on galaxies making the structures in 
the disk plane difficult to detect. And very few galaxies 
have the pc-scale quality of kinematic data that we have 
here. However, where this data exists there is evidence 
for nuclear structures similar to those presented here. 

The Milky Way hosts a nuclear star cluster centere d 
on Sgr A*, first noted by iBecklin fe Neugebauerl (|1968l) . 
The luminosity profile of the inner ^100 pc of the Milky 
Way is distinc t from the und erlying bulge and disk pro- 
files (e.g. ISerabvn fc Morris 1996). This component is 
clear ly flattened along the galactic plane (axial ratio of 
0.45; ICatchpple et~ al. 1990), b ut less so in the central 30 
pc, which led lLaunhardt et a l. (2002) to propose the exis- 
tence of two components, a 3xl0 7 Msol NSC dominating 
the central 30 pc (with r e ff ~ 10 pc) and a more massive 
nuclear stellar disk extending out to 100-200 pc. Numer- 
ous young OH/IR stars are found, primarily within the 
central 40 pc and these ha ve a flattened distribut ion with 
rotation of >100 kins" 1 (Lindqvist et al. 1992). Based 
on integrated K-band spectra, McGinn et al. (1989) find 
that the central dispersion is ^125 kms -1 and drops 
quickly with increasing radius, while the rotation in- 
creases outwards equalling the rotation at a radius of 
^3 pc, the edge of their study. Compared to the 
NGC 4244 NSC, the Milky Way appears quite similar, 
but scaled up by ~10x in mass and ~3x in size, ro- 
tation, and dispersion as would be expected for a more 
massive galaxy. The very center (<1 pc) of the MW NSC 
has been intensively studied and shows complex dynam- 
ics with striking differences between older stars and very 
young (<10 Myr) stars. The latter show ordered rotation 
(and cou nter rotation!) an d provide ev idence for pro- 
cesses 2 JG cnzcl ct al. 2003) and 3 (e.g. Maillar d et al] 
12004 ILu et al.ll2005h described above. 

There is also evidence for flattening and rotation in the 
M33 NSC. This cluster is elongated along th e major axis 
of th e galaxy with an axial rat io of ~0.85 (|Lauer et alj 
H99l iMatthews fe Wood! l200l . At a radius of ~1" 
(4 pc), th e rotation is ~8 km s -1 and dispersion is 
27kms _1 (jGebhardt et al.ll200lh . Even considering that 
the orientation is likely not edge-on, it appears that the 
rotation does not dominate in the M33 NSC to the same 
extent as in NGC 4244. 
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In conclusion, we have shown above that the nuclear 
star cluster in NGC 4244 is rotating. The rotation is seen 
in both the young disk and older spheroid components, 
suggesting that accretion of gas or star clusters from the 
disk is the dominant mechanism in the formation of this 
NSC. While there is some evidence that other NSCs in 
spiral galaxies have similar structures (and thus forma- 
tion histories), a better understanding of NSC formation 
across Hubble type will require a detailed study resolv- 
ing the stellar populations and kinematics of numerous 
nearby NSCs. 
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